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A simple, practical model for computing the equilibrium thermodynamics and structure of jellium 
by classical strong coupling methods is proposed. An effective pair potential and coupling constant 
are introduced, incorporating the ideal gas, low density, and weak coupling quantum limits. The 
resulting parameter free, analytic model is illustrated by the calculation of the pair correlation 
function over a wide range of temperatures and densities via strong coupling classical liquid state 
theory. The results compare favorably with the first finite temperature restricted path integral 
Monte Carlo simulations reported recently. 
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The limitations of many body theories for strongly cou- 
pled quantum systems at finite temperatures have led to 
attempts to adapt corresponding methods known to be 
effective for classical systems [1]. Among these are the 
classical molecular dynamics (MD) simulation method, 
classical Monte Carlo integration, and liquid state the- 
ory [5] , modified with effective pair potentials that incor- 
porate essential quantum effects such as diffraction and 
degeneracy. Early approaches were based on a classical 
form for the two particle density matrix in coordinate 
representation to identify the effective pair potential in- 
corporating diffraction effects [21 E]. Exchange effects 
were incorporated in a similar way using the pair corre- 
lation function for an ideal gas [5l|6]. More recently, effec- 
tive classical systems have been defined with an effective 
temperature as well as pair potential [SI E]. A formal- 
ism for construction of a classical system with thermo- 
dynamics and structure corresponding to a given quan- 
tum system is described in reference [8j. A system of 
particular interest exhibiting strong Coulomb coupling 
effects, both classical and quantum, is the electron gas 
(referred to classically as the one component plasma or 
quantum mechanically as jellium). In the classical limit 
its thermodynamics is completely characterized by the 
Coulomb pair potential and a dimensionless coupling 
constant F = Pq^/rQ. Here (3 = l/ksT is the inverse 
temperature, q is the particle charge, and tq is the av- 
erage distance between particles defined in terms of the 
density (see below). 

The objective here is to provide an effective classi- 
cal system representing the thermodynamics of the real 
quantum jellium, using an effective pair potential and 
an effective coupling constant. Simple analytic expres- 
sions are given, based on the more complete but complex 
results of reference [S] . Application of this model is illus- 
trated using the hypernetted chain (HNC) integral equa- 
tion of classical liquid state theory to calculate the pair 
correlation function. Comparison of these calculations 
with the first finite temperature restricted path integral 
Monte Carlo (PIMC) simulation results reported recently 



[lOj show good agreement over a wide range of densities 
and temperatures. 

Jellium has a broader current interest than its histor- 
ical role as a test bed for quantum many body theo- 
ries. Its thermodynamic properties, particularly the den- 
sity dependence of the free energy, provide the basis for 
local density approximations (LDA) in electron density 
functional theory (DFT) [11]. Fits to zero temperature 
PIMC simulations have been the basis for virtually all 
LDA DFT for the past thirty years. Until now there has 
been no corresponding basis for an LDA at finite tem- 
peratures from either theory or simulation. Such con- 
ditions of solid densities at temperatures comparable to 
the Fermi temperature are of central interest to the new 
studies of "warm, dense, matter" [T^ . 

The system of interest is a collection of N charges with 
Coulomb pair interactions </)(r) in a uniform neutralizing 
background, at equilibrium |l3j . The thermodynamic 
variables are the temperature and density, T = l/ksP 
and n. A corresponding effective classical system is con- 
sidered with pair interactions 4)c{t) in a uniform neu- 
tralizing background, at equilibrium with temperature 
and density Tc and Uc- The correspondence of the clas- 
sical and quantum systems is established by defining 
(j)c{r),Tc,nc as functions or functionals of (j3{r),T,n in 
such a way as to assure the equivalence of selected equi- 
librium properties. Three such conditions are chosen [5]. 
The first two are equivalence of the densities and pair 
correlation functions 

nc = n, gc{r,l3c,nc\ (l)c) = g{r,l3,n\ (j)). (1) 

The remaining condition fixing Tc is replaced by a corre- 
sponding condition for an effective coupling constant, as 
discussed below. 

To be useful, the condition equating pair correlation 
functions must be invertible, ipcir) — g~^{r, l3c,nc \ g), 
which entails solution to the classical many-body prob- 
lem, (this inversion does not need to be unique; see fi- 
nal comments below). In the special case of the ideal 
gas limit the result is known as the Pauli potential. 
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denoted 4>^c'^ (r). Even in this case the inversion can- 
not be accomphshed exactly, but good approximations 
are known [5J |2] • The relevant dimensionless thermody- 
namic parameters for the quantum system are the tem- 
perature relative to the Fermi temperature t = l/Pep 
and rg ~ ra/aB, the mean distance between particles 
To (defined by AirnrQ/S — 1) relative to the Bohr ra- 
dius qb- Hence the dimensionless effective potential 
(j)*{x, t, Tg) = I5c4>c{f, n I 0) is written in the form 



<^:(x, t, r,) = (x, t) + A* (x, t, r,) , 



(2) 



where x = rjrQ. It has been recognized that the Pauli 
potential depends only on t. 

Two exact limits for A* (x, i, r^) are important for the 
discussion here. The first is the weak coupling limit 



(3) 



where Cc(s,t, r^) is the direct correlation function. It 
is related to the pair correlation function gJ^x^t^Tg) by 
the exact Ornstein-Zernicke equation [5]. Using the cor- 
respondence conditions ([T]) the Ornstein-Zernicke equa- 
tion defines the direct correlation function in terms of the 
quantum pair correlation function 



Cc[x) = g{x) - 1 
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dx'c(|x-x'|)[.9(a;')-l]. (4) 



A sufficient condition for weak coupling is large x, for 
which the behavior of Cc{x,t,rs) is determined from the 
perfect screening sum rule for g{x,t,rs) [H], giving 



A* {x,t,rs) it,rs)x-\ 



(5) 



This is the same form as for the classical one component 
plasma, except with the classical Coulomb coupling con- 
stant r = /Sq'^/ro = 4 (2/37r2)^/V,/3i replaced by the 
effective coupling constant 



Tc {t,r,) 



Phuip coth {(3h(jjp/2) 



(6) 
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Here ujp = (47mg^/m) is the plasma frequency (or, 

equivalently, fShojp = 4 {2V3tt-'^)^^^ rV^ /3t). At low 
temperatures and fixed density T becomes divergent 
whereas the effective coupling constant remains finite 

1/2 

Tc (0, Tg) — 1. 155 rs . At high temperatures Tc (t, Vg) — ?> 
r ~ 0. 543rg/t. 

The second exact limit is that for low density and 
weak coupling. The condition of low density means that 
g{x, , t, Vg) is determined by the two electron Slater sum. 
The weak coupling A* (x,t,rs) in that case is known as 
the Kelbg potential [3ll3tfT4] 

A* {x,t,rg) Tx"^ (l - exp(- (ax)^) 

+y/TT{ax)erfc{ax)) , (7) 
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where a — (rg/T) . This weak coupling result at low 
density can be improved by imposing the exact behavior 
of the two particle Slater sum at a; = 0, to include some 
strong coupling effects [31 IH [151 [E] • The modified form 
is 



exp((aa;) ) 



/tt — erfc{^ax) , (8) 
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with 



7 (rr^) = 



(TrTrg) 



1/2 



lnS'(rr, 



(9) 



where S{Trg) is the two electron relative coordinate 
Slater sum at a; = 



SiTrg) = -4(^rr,)^/' / dye-y 



y 



o7r(rr. 



-,1/2 



/y 



(10) 

The proposal here is to further extend this Kelbg form 
to apply broadly across a wide range of values t,rg by 
imposing the exact asymptotic limit ([5|. This is accom- 
plished by replacing F with the effective coupling con- 
stant Fc given by ([6]). The approximate effective pair 
potential is thus 



0:(x, t, rg) ~ (x, t) + A*j, (x, F„ rg) . 



(11) 



^ {x,Tc,rg) is an analytic, parameter free form 



Since A^^ (x,F 
this potential is suitable for practical applications in clas- 
sical many-body theory, classical Monte Carlo calcula- 
tions, and molecular dynamics simulations. 

To illustrate the utility of this model potential the pair 
correlation function g{x,t,rs) for jellium is calculated 
here using the classical liquid state HNC integral equa- 
tion [2 . The first step is a determination of ^c*^-** (x, t) for 
the ideal Fermi gas. Since the pair correlation function 
g^^'' (x, t, rg) is known exactly, the HNC equations can be 
inverted to determine (j)^^^* {x,t). These equations are 
solved numerically using the method of Ng [T^]. Next, 
with (/)c^^*(x,t) known the pair correlation function for 
jcUium can be determined from the HNC equations us- 
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Very recently restricted path integral Monte Carlo 
(PIMC) simulations have been reported for the pair cor- 
relation function spanning conditions ranging from ex- 
treme quantum to semi-classical. These results provide 
important benchmarks for existing quantum many-body 
methods, as well as the approach proposed here. Con- 
sider first the relatively high temperature t = 8. Figure 
(la) shows good agreement with PIMC for all densities, 
1 < rg < 40. Figure (lb) shows the same comparison 
for i — 1. Again the agreement is good, except at the 
extreme condition rg — 40. In this case a strong corre- 
lation peak has formed that is badly under estimated by 
the theory, although its location is adequately described. 
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FIG. 1: Pair correlation function g{r*) at (a) t = 8 and (b) t = 1 for = 1, 6, 10, 40. Also shown are the results of PIMC. 



Generally, it is found for t < I the theory is quite good 
for 1 < Ts < 10. This is illustrated in Figures (2a) and 
(2b) a.t t — 0.5 and 0.0625 (the latter is essentially the 
same as i = 0, as confirmed by a comparison with diffu- 
sion Monte Carlo simulations at t = [20]). Some trends 
are evident even from this limited data. For example, 
the temperature dependence for t < 0.5 is quite weak for 
Ts > 1. However, for rg — 1 a, significant temperature de- 
pendence is seen for 0.5 < i < 8. In summary, the model 



potential (111 provides a practical form for the analysis 



of jellium using classical methods under conditions that 
are difficult to access by existing quantum methods (e.g., 
Ts > 1 and t < 10). 

The thermodynamic properties of jellium can be cal- 
culated from the pair correlation function. For example, 
the pressure can be obtained from a coupling constant 
integration. Let p{t,rs,q) be the exact quantum pres- 
sure and g (r, i, Ts, q) the exact quantum pair correlation 
function where now the dependence on the charge q has 
been made explicit. Then 

p{t,rs,q) = p(i,rs,0) + 87r / dyy 

Jo 

/"OO 

X / drr^^{r){g{r,t,rs,y)-l). (12) 
Jo 

Here (j){r) is the Coulomb pair potential. Therefore, ap- 
proximating g (r,t,rs,y) by the corresponding classical 
result obtained using the model potential ( 11 ) determines 
the pressure for arbitrary t,rs. A more direct approach 
would be classical Monte Carlo integration of the Gibbs 
distribution for the free energy 



-i:(<A<'"*(^M,*)+AK(^>.,r„r,)) 
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(13) 



with Xij = |xi — Xjl . 

As noted in the introduction, the idea of an effective 
classical pair potential with an effective classical temper- 
ature was already introduced more than ten years ago by 
Perrot and Dharma-wardana [6]. Instead of the Kelbg 
potential they chose the Deutsch potential [17 , originally 
introduced as a simpler representation of the Kelbg po- 
tential. The PDW effective classical potential is similar 
to (111, but with {x,Vc,rs) replaced by 



A 



PDW 



{x,TpDw,rs) ^TpDwx ^ (1 - exp(-6x)) (14) 
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Here, b = {'nrg/TpDw) and the effective coupling con- 
stant is 



PDW 




-1/2 



(15) 



This follows from their phenomcnological form for the 
classical temperature interpolating between the real tem- 
perature T and a finite value Tq at T = 0. The sin- 
gle parameter Tq/T is determined by fitting the classical 
correlation energy calculated with this potential to the 
quantum exchange/correlation energy determined from 
PIMC at T=0. It is given as an explicit fitting func- 
tion of Ts in reference [B]. Although the dependence of 
TpDW on t,rs is quite different from that derived here, 
and the shape of the resulting effective pair potential can 
be quite different, nevertheless the HNC pair correlation 
function calculated from the PDW potential has a sim- 
ilar accuracy to that reported here. This indicates that 
an effective pair potential has no inherent physical inter- 
pretation, but rather is a non-unique tool for generating 
physical properties of interest through classical many- 
body methods. Here that potential has been constructed 
by imposing three exact constraints: the ideal gas limit. 
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(a) t=0.5 




X 



FIG. 2: Pair correlation function g{r*) at (a) t — 0.5 and (b) t 

low density limit, and large distance limit. Consequently 
no fitting parameters are required. The result provides 
theoretical support for the ideas of reference and pro- 
vides insight into the relevant physical mechanisms. For 
example, the exact screening sum rule that determines 
the form of Tc here appears to incorporate quantum ef- 
fects as significant as those of Tp]j\Y imposed by empir- 
ical T — exchange/correlation energy data. 
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